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Abstract: We introduce an alternative parametrisation for the scale dependence 
of the non-linearity parameter /nl in quasi-local models of non-Gaussianity. Our 
parametrisation remains valid when /nl changes sign, unlike the commonly adopted 
power law ansatz /nl(&) ck k nf nl. We motivate our alternative parametrisation by 
appealing to the self-interacting curvaton scenario, and as an application, we apply it to 
the CMB power asymmetry. Explaining the power asymmetry requires a strongly scale 
dependent non-Gaussianity. We show that regimes of model parameter space where /nl is 
strongly scale dependent are typically associated with a large <?nl and quadrupolar power 
asymmetry, which can be ruled out by existing observational constraints. 
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1. Introduction 

The paradigm of inflation [1, 2, 3, 4, 5] is the simplest framework which explains the origin 
of the primordial cosmological fluctuations [6, 7, 8, 9, 10]. The detailed information which 
is encoded in these fluctuations over a range of cosmological scales offers the opportunity 
to test the plethora of inflationary models which have been proposed in the literature. 
In particular, the local ansatz for non-Gaussianity has become a popular description for 
modelling deviations from a purely Gaussian distribution of primordial perturbations. 

For the simplest model where the primordial perturbation is a quadratic local function 
of a single Gaussian field, non-Gaussianity is usually parametrised in terms of a single 
constant parameter, /nl, which measures the amplitude of the bispectrum normalised to 
the square of the power spectrum of primordial curvature fluctuations [11]. In this case /nl 
is scale-independent by construction. In general however /nl is not constant over all scales, 
as expected from theoretical considerations [12, 13, 14, 15, 16, 17] (see [18] for considerations 
of scale dependent equilateral non-Gaussianity). In this paper we exclusively consider 
non-Gaussianity which reduces to the local shape in the limit of scale invariance. Scale- 
dependence can arise due to multiple fields contributing to the curvature perturbation, a 
curved field space metric [19, 20], or due to non-linear evolution of modes after they leave 
the Hubble radius during inflation. This scale dependence is typically parametrised as a 
power law, 


/nl(&) — /NL(^’/^o) n/NL 


(1.1) 




where /^ L is the amplitude of /nl measured at some pivot scale ko and n/ NL is a constant 
(see, e.g., [21, 13, 22] for details). Re-writing this as 


_ dln|/ NL | 
n/NL_ d Ink 


( 1 . 2 ) 


one can clearly see that n/ NL diverges whenever /nl changes sign. The self-interacting 
curvaton model [23, 24] and the dynamics of perturbative reheating in the axion-like model 
studied in Ref. [25] are two examples where /nl can change sign over the observable range 
of scales, 10 3 Kpc < k~ 1 < 10 4 Mpc. In such examples the scale-dependence of /nl cannot 
be described by a simple power law over scales which include /nl close to zero. We define 
/nl as the ratio of the bispectrum to the square of the power spectrum in the usual way, 
which for an equilateral configuration reduces to 


/nl(^) 


5 B^(k,k,k) 
18 P*(h) 


(1.3) 


In this paper we introduce an alternative parametrisation for the scale dependence of 
/nl (k): 

fNh(k) = /nl + 7/nl l n (/[[) ^ -’ (l-^) 

where q/ NL is a constant, and /^ L is the value of /nl measured at the pivot scale ko for 
an equilateral configuration. Unlike the conventional parametrisation /nl oc k nf nl . ours 
remains valid when /nl vanishes, and is theoretically well-motivated as we demonstrate by 
appealing to the self-interacting curvaton model as an explicit example where /nl(&) can 
exhibit a sign change over an observable range of scales. We also show that even in cases 
with |7/ nl | 1, it is often a good approximation to take this parameter as a constant over 

the observable range of scales. 


As an important application of our results, we study the observed dipolar power asym¬ 
metry of the primordial power spectrum [26]. On large scales this is observed to be a large 
effect, but the asymmetry must be at least an order of magnitude smaller on small scales, 
where the effect is not observed but tightly constrained [27]. If not a statistical fluctuation, 
the asymmetry may arise through the non-Gaussian coupling between a large amplitude 
super horizon mode and the modes where the asymmetry is observed. The amplitude of the 
non-Gaussianity must also vary by an order of magnitude in order to explain the reduction 
of this effect on small scales. Such a large variation cannot be accurately described by the 
conventional parametrisation /nl(^) = /NL(^Ao) n/NL with constant n/ NL , which is only 
valid under the assumption |n/ NL ln(k/ko)\ <C 1, however we explicitly show that it can be 
accurately described by Eq. (1.4). 

We also show that regions of model parameter space where /nl varies sufficiently 
quickly with scale and is large on the scales where the asymmetry is large typically give 
rise to a large c/nl and quadrupolar modulation of the power spectrum. Unlike /nl and 
the quadrupolar asymmetry, these two additional effects are close to scale invariant and for 
typical parameter values too large to be consistent with existing observational constraints. 
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In the next section we introduce our parametrisation of strongly scale dependent /nl 
and apply it to the self-interacting curvaton scenario in Sec. 2 . 1 , deriving new analytic 
solutions in Sec. 2.2 which are valid in the limit that the quadratic term of the curvaton 
potential dominates. Unlike previous analytic solutions, ours takes into account the evolu¬ 
tion of the curvaton during inflation and reheating, and we show that this evolution has a 
large effect which may give rise to qualitatively new behaviour of the non-Gaussianity. We 
apply our formalism to the power asymmetry in Sec. 3 and study the difficulty in keeping 
( 7 nl and the quadrupolar power asymmetry small in Sec. 3.1. Finally we conclude in Sec. 4. 
Throughout this paper we work in units where h = c = 1, and we set the reduced Planck 
mass M p = 1. 

2. Strongly scale dependent non—Gaussianity 

During inflation, it is normally accurate to parametrise the power spectrum by a power 
law, both from a theoretical and observational perspective. For both single and multi-field 
inflation, it is normally accurate to treat the spectral index as a constant over a much 
larger range of scales than we can observe (unless the running is large). The spectral index 
is related to the slow-roll parameters, which are normally required to be small in order 
for inflation to last many efoldings. A spectral index of local /nl, n/ NL , was introduced 
in [28, 12, 21, 13, 22] in analogy with the power spectrum and its value is also proportional to 
the slow-roll parameters. However, in contrast to the power spectrum, the scale dependence 
of /nl can be greatly enhanced by numerical factors much greater than unity [24] . In such 
cases, it is not accurate to treat n/ NL as a constant, but it is precisely such cases which 
are the most observationally interesting because the current tight constraints on /nl mean 
we can only hope to measure both /nl and its scale dependence if the scale dependence is 
much larger than the observed scale dependence of the power spectrum, for forecasts see 
e.g. [21, 29, 30, 31, 32], 

Motivated by the observation that regions of large ny NL tend to correspond to regimes 
where /nl is small, we introduce the new parametrisation Eq. (1.4) which remains valid 
as /nl crosses zero. Provided that /nl grows to large values within the scales that we can 
probe, this case is observationally interesting. As an additional motivation, the observed 
power asymmetry which we discuss in Sec. 3 requires a strongly scale dependent /nl, with 
n/ NL ~ —0.5. For most models, treating n/ NL as constant will only be valid for about 2 
e-foldings when its magnitude is so large, which is a much smaller range of scales than over 
which the observation is made. Our parametrisation can be easily related to the standard 
one 

7/nl = = /nl( fc)n /NL (fc), (2.1) 

but our definition remains valid in the limit /^ L = 0. Furthermore we will see that 
our parametrisation remains valid (for constant 7 / NL ) over a large range of scales in the 
observationally interesting regime, which includes the requirement that |t/ nl I 1 - In 
contrast, in such regimes where /nl is strongly scale dependent, it is clear from (2.1) that 
rif Nh must be equally strongly scale dependent in order for 7 y NL to be a constant. 
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2.1 The self-interacting curvaton scenario 

In this section we motivate the parametrisation Eq. (1.4) by studying the scale dependence 
of /nl in the self-interacting curvaton scenario. The model is described by the potential 

a) = V(<f>) + l -m 2 y + \a n , (2.2) 


where a is the curvaton field, and inflation is driven by V ((/>)■ We will consider n > 4, 
with n even. The phenomenology of this model has been extensively studied in previous 
work, see e.g. [33, 23, 34, 35, 14], where it was found that observable quantities (such 
as /nl) can depart substantially from the standard quadratic curvaton prediction due to 
self interactions. It is useful to define the following quantity 


$k — s(tk) 


2A a n -\t k ) 


m- 


(2.3) 


which gauges the strength of the curvaton self interaction relative to the mass term in the 
potential. Throughout this paper, we denote the time at which the mode k left the Hubble 
radius by tk■ For the sake of brevity we will often write Xk = X(t k ) for scale dependent 
quantities which are to be evaluated at the epoch of Hubble exit. 

After inflation has ended (f> begins to oscillate about the minimum of its potential, 
decaying to radiation as it does so, which comes to dominate the energy density of the 
universe. We employ the sudden decay approximation to model the decay of the curvaton 
field, which assumes that the curvaton decays instantly into radiation when H = T a . While 
the sudden decay approximation does not properly treat the gradual transition from an 
oscillating scalar field to a radiation bath, the behaviour far from this transition period, 
and the resulting ratios of energy densities are well captured by the approximation. As has 
been shown in previous work [36, 25] this is enough to quite accurately reproduce the effects 
of the decay process on primordial observables. If the curvaton is coupled directly to the 
inflaton or to other scalar Helds, then it can decay non perturbatively through parametric 
resonance [37, 38, 39, 40]. Since in the potential Eq. (2.2) the fields are coupled only through 
gravity, we do not consider this possibility here. Sufficiently deep into the phase of coherent 
oscillations of the curvaton field, but before it has decayed into radiation, the curvaton 
energy density is very well described by the quadratic contribution, p a = rn^a 2 sc /2, where 
CT OS c = a osc (a(t k )) is the amplitude of the curvaton oscillations about the quadratic part 
of the potential. We assume that the curvature perturbation ( is generated solely by 
fluctuations of the curvaton field and we neglect the inflaton fluctuations. Appealing to 
the 5N formalism [41, 42, 43], the curvature perturbation can be written at third order as 
(see e.g. [36]) 


Ck = N'(t k )5a k (t k ) + }-N"(t k )(5a -k dcr) k (4) + ^N'"(t k )(6a -k 5a * 5a) k (t k ) + ... 
2 fa 


O (7 ns r. O 


+ 


T dec 


^osc ( 

^osc 

jn 

osc 


// 

^OSC | 

( c4cV 

°osc 

V ^osc / 


(7, 


+ 3 


// / 

• ^OSC^OSC 


(7n 


(7Z 


(5a * da ) k (tfc) 
(5a -kSa-k 5a) k (tk ) + ... , 


(2.4) 
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where 


3 Pa 


(2.5) 


2>p a + 4p 7 dec 


Subscript ‘dec’ denotes the time of curvaton decay, and (') = d ak . The convolutions are 
defined by (5a * 6a)^(tk) = (27r)~ 3 / dq(kr q (tfc)<Icxk-q(£fc)- N(tk) denotes the number of e- 


foldings from an initial spatially-flat hypersurface at tk, to a final uniform density surface 
after the decay of the curvaton. The non-linearity parameters /nl(&) and 5nl(^) are given 
by (see e.g. [36]) 



( 2 . 6 ) 


Scale dependence is manifest since < 7 osc is a function of a^. If the curvaton potential were 
exactly quadratic, the equation of motion for a would be linear, and / osc = 1 , g osc = 0 . 
When the potential deviates from quadratic, derivatives of a osc contain information about 
the self interaction. As has been previously observed [14, 23], even a small correction to 
the quadratic potential can have an important effect on the level of non-Gaussianity while 
leaving the Gaussian perturbation essentially unchanged. 

In the limit r^ ec ^ 1) we see from Eq. (2.6) that /nl vanishes whenever f osc = 
4rdec/3 <C 1. Working in this limit and using the standard parametrisation for fNh(k), 
Eq. (1.2), the authors of Ref. [24] studied the scale dependence of /nl ( a nd 5nl)- For 
n > 6 , /nl wa s found to be a highly oscillatory function of 07 ., resulting in a series of 
‘spikes’ in n/ NL as /nl changes sign. At these points, the standard parametrisation cannot 
be used. 

2.2 Scale dependent /nl in the limit < 1 

Analytic solutions in the limit of weak self-interactions have been obtained in previous 
studies [23] however, these solutions are restricted to the regime where the curvaton remains 
frozen until the inflaton has decayed into radiation. The solution that we present below 
accounts for evolution of the curvaton field during inflation. One motivation for including 
the effects of curvaton dynamics during inflation is that these dynamics may be relevant 
for the CMB power asymmetry. Because we are interested in the parameter space with 
large non-Gaussianity, i.e. r& ec < 1, we assume throughout this paper that the curvaton 
does not effect the Hubble parameter. We have numerically checked that this is a good 
approximation. 

In the limit Sk -C 1 the potential Eq. (2.2) is dominated by the mass term, and the effect 
of the curvaton self interaction is to weakly deform the potential away from pure quadratic 
form. As the curvaton rolls down its potential, the magnitude of the self-interaction term 
relative to the mass term quickly decreases, until a point is reached where the potential is 
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almost exactly quadratic. In this regime, the equation of motion for a is linear, 

a + 3H& + rn 2 a a « 0 , for t > t q . (2-7) 

The general solution can be expressed in the form a(t) = a osc [aY\ (t) + 6 V 2 (t)] where <t osc = 
a(t = t q ), and t q is defined to be the time from which the linear equation above becomes 
valid. If we make the additional assumption that a does not contribute to the Hubble rate, 
Yi, 2 (t) are functions of time which depend only on the inflationary background. Owing to 
the non-linearity of the curvaton equation of motion for t < t q , a osc will be a non-linear 
function of a k . To calculate this functional dependence we assume that a remains slowly 
rolling until Eq. (2.7) becomes a faithful description of the curvaton dynamics (note that if 
the curvaton is sufficiently light, it will remain in slow-roll even after inflation ends). This 
enables us to write a osc = /3asR(t = t q ), for some constant (5, where <tsr is the solution 
obtained from the curvaton slow-roll equation of motion: 

3H& + m^cr + nA(X ri_1 « 0, for t<t q . (2.8) 


This equation can be directly integrated, to yield the solution at time t 
inflationary background: 

1 

e -fj tT (n-2)l(t q ) n -2 

1 + f Sfe [l - e -V*(.n~2)Z(t q )] 


&SR{t q ) = (T k 


t q for a generic 


(2.9) 


where 

Sl'wr 

and the parameter fj a = is related to the curvaton slow-roll parameter 


rjaitk) 


W k 

3Hl = ?)<T 


1 H— n(n 
2 v 



( 2 . 10 ) 


( 2 . 11 ) 


The curvaton does not contribute to the integral 1(t q ), since we have neglected its back- 
reaction to gravity. To a good approximation, the curvaton begins to oscillate when the 
Hubble rate becomes equal to the curvaton mass. An appropriate choice for t q is there¬ 
fore H 2 (t q ) ~ m The integral X(t q ) is dominated by its upper limit, and so we have 
approximately I(t q ) ~ ~ 1 ff] a . Inserting this result into Eq. (2.9) and making the 

simplification e~( n ~ 2 > ~ 0 for n > 4, we have: 




= /3ctsr ~ /3e ° w a k (l + ^s k ^j 


( 2 . 12 ) 


Taking derivatives with respect to a k and inserting into Eq. (2.6) we obtain the remarkably 
simple result for /nl in the limit r^ec ^ 1 : 


fNh(k) 


5 

d^dec 




5 

3 ' 


(2.13) 


This expression relates /nl to the curvaton self-interaction strength s k , which was defined 
in Eq. (2.3). The scale dependence of /nl is manifest since s k is a function of k. We 
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Figure 1: The asymptotic value of /nl (the final value which should be compared to observations) 
as a function of the self-interaction strength Sk- We have chosen an arbitrary value r<j e c = 0.05. 
Solid lines correspond to a numerical calculation, and dashed lines to the analytic result Eq. (2.13). 
For the numerical calculation, we took V(<f) = rri^cj) 2 / 2, with </>& = 16, 'ij a = 2.3 x HR 2 . For this 
and all other figures, we choose T^ = H 0sc , where H osc is the Hubble rate at the time the inflation 
starts to oscillate. The intersection of the horizontal dashed line with the analytic solution indicate 
the points, Sk = Sq, where /nl = —5/3. These values are given by Eq. (2.16). 


compute this scale dependence explicitly in what follows. When Sk = 0, we recover the 
standard quadratic curvaton result, /nl = 5/(4rd ec ) — 5/3, which is independent of k. 
By virtue of the approximation e-Vcr(n-2)i(t q ) ^ our resu p, f or / NL hag no dependence 
on fj a . Had we computed the integral I(t q ) exactly (which would require specifying the 
inflationary potential V ((/)), /nl would develop a weak dependence on fj a . In the limit 
S|;<Cl this dependence is negligible, and so to a very good approximation /nl is a function 
of Sk only. To the best of our knowledge, this is the first time an analytic solution for /nl 
- which is valid for arbitrary fj a - has been constructed for the self-interacting curvaton 
model. An analytic solution for /nl in the limit Sk <C 1 was obtained in Ref. [23], and 
a solution for n = 4 valid for all Sk was presented in Ref. [44]. Both of these solutions 
however are restricted to the regime where the curvaton remains frozen until the inflaton 
has decayed into radiation. To consistently realise this scenario within a typical inflationary 
background requires fj a <C HR 7 . As we shall see shortly, the scale dependence of /nl is 
directly proportional to fj a , and so to generate an interesting (large) scale dependence 
requires fj a 10 -7 . 

In Fig. 1 we compare Eq. (2.13) as a function of the self-interaction strength Sk, against 
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Figure 2: The asymptotic value of /nl as a function of the self interaction strength s*,, for n = 4, 
for various values of fj a . For large values of fj a , / NL can cross zero. For fj^ = 10 -8 (black line) 
the curvaton remains frozen until the inflaton has decayed, and we recover the results of Ref. [24]. 
For Sk 1, /nl is independent of fj a , as revealed by the analytic solution Eq. (2.13). All lines 
are computed numerically, with V(4>) = rrv^<p 2 /2, <pk = 16. We have chosen an arbitrary value 
r dec = 0.05. 

a full numerical calculation, for various values of n. Our numerical technique (which is de¬ 
scribed in detail in Ref. [25]) computes the derivatives appearing in Eq. (2.6) for /nl 
exactly, accounting for gravitational backreaction of the curvaton. For illustration pur¬ 
poses, we choose a quadratic potential for the inflaton, V(4 >) = nri^cj) 2 / 2, however we have 
confirmed that the results for /nl are largely insensitive to this choice. In our numerical 
calculations, the inflaton transfers its energy to radiation through a friction term T^cj), (i.e., 
<j) + (3 H + r 0)0 + = 0 ) which we introduce by hand into the equation of motion once 

</ is coherently oscillating. Throughout this paper, we set T^ = H osc , where H osc is the 
value of the Hubble parameter at the onset of coherent oscillations of the inflaton field. 
With this value of T^, the inflaton decays to radiation almost instantly. We have checked 
numerically that the results are largely insensitive to this choice. 

It is possible for /nl to change sign for n = 4 when the curvaton is allowed to roll 
during inflation, as clearly seen in Fig. 1. This behaviour was not observed in earlier 
studies of the self-interacting curvaton model with n = 4 [23, 24] which were restricted 
to the regime where the curvaton remains frozen throughout inflation. We do however 
numerically recover the results of [24] for n = 4 as fj a is steadily decreased until the point 
where the curvaton remains frozen throughout inflation. This is illustrated in Fig. 2. 
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As can be seen from Fig. 1, /nl changes sign as the self-interaction strength is varied. 
The standard parametrisation for /nl, Eq. (1.2), is not valid if this sign changes occurs 
over the observable range of scales. The alternative parametrisation /niX^) oc In A: for 
/nl = i- e - Eq. (1.4) does however remain valid. We now make use of Eq. (2.13) for /nl 
to motivate this alternative parametrisation. 


Our objective is to expand Eq. (2.13) for /nl about some pivot scale k = ko, which we 
choose to be the scale at which /nl = —5/3. Fluctuations on observable scales 10 3 Kpc < 
A: -1 1$ 10 4 Mpc today correspond roughly to A N ~ 8 e-foldings over which these modes 
left the Hubble radius during inflation. Over this range of N , the inflationary background 
is a slowly varying function of time. Hence, the integral X(t) given by Eq. (2.10) may be 
expanded about N ~ In (k/ko) = 0. To first order in slow-roll we obtain: 


Xsr 


In (k/ko) 1 + e^ln (k/k 0 ) 


(2.14) 


where e<^ 0 = —Hq/Hq, and a subscript ‘0’ denotes evaluation at the pivot scale ko. Substi¬ 
tuting the above result for Zsr. into the slow-roll solution for the curvaton, Eq. (2.9), and 
using Eq. (2.3) for Sk yields: 


&k So 


e ~Vcr 0 (n—2 )X S r 

1 + [1 - e-^oO-^sR] ’ 


where sq is the value of s/- at the pivot scale: 


2 

s ° = “7 - 

n(n — 1) 


(2.15) 


(2.16) 


which implies r] ao = 2r/ CTo , see Eq. (2.11). Finally, substituting Eq. (2.15) into Eq. (2.13) 
for /nl and expanding the result to second order in In (k/ko) we obtain: 


/nl(*0 = _| + ^ a? l r7 In (k/k 0 ) 
3 8r dec 


1 -%(^+ 1 ) (va 0 - j— v 

4 \ a n [n H - 1) 


e^o In (k/ko) 


, (2.17) 


where a n = (n — 2)/(n — 1), and we have used r] ao = 2fj (T0 . The scale dependent part 
of /nl vanishes for n = 2 (where a n = 0), which is simply the non-interacting limit of 
the model. The constant term and the term linear in Yn.{k/ko) correspond exactly to the 
parametrisation of /nl we introduced in Eq. (1.4), with 

f o _ 5 _ 5na n 

/nl o ’ 7/nl Q„ • (2-18) 

4 oT d ec 

In Fig. 3 we compare Eq. (2.17) against a fully numerical calculation for two different values 
of n and r] aQ . For r/ CT0 > 10 -2 , the term quadratic in In (k/ko) becomes important, and the 
scale dependence becomes non-linear in lnfc. 
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Figure 3: The scale dependence of /nl- Solid lines correspond to a numerical calculation, and 
dashed lines to our analytic result, Eq. (1.4). The left panel corresponds to n = 6 and the right 
panel to n = 8. Note the different scales on the y-axis. For r/ ao > 0( 10~ 2 ), /nl becomes a non¬ 
linear function of ln(fc/fc 0 ). For the numerical calculation, we choose V{<j>) = 2, <j>k = 16. 

We have chosen an arbitrary value rdec = 0.05. 

3. The CMB Power Asymmetry 

Both WMAP and Planck have observed a dipole asymmetry in the power spectrum on 
large scales at about 3a significance. This is parametrised by 


P( = P iS o(l + 2Tfi.p) , (3.1) 

where Pj so is the isotropic power spectrum, n is a direction on the sky of ( and p is 
the direction in which the asymmetry is maximised. The current best measurement is 
A = 0.06 — 0.07 [45] which is based on measurements on large CMB scales, 1 = 2 — 64, 
while the smaller scales show no evidence for an asymmetry, requiring A < 0.0045 for 
l = 601 — 2048 [46, 47, 48]. Explaining the asymmetry from an inflationary model requires 
both a large amplitude mode with wavelength larger than the horizon scale today and non- 
Gaussianity to correlate this very large scale mode with the CMB scale modes which are 
asymmetric [49, 50, 51]. The super horizon scale mode must originally be an isocurvature 
perturbation and cannot be generated from single-field inflation [52, 53]. It’s existence is 
typically just postulated, but might be related to the beginning of inflation [54] or global 
curvature [55]. The asymmetry could also be associated to a domain wall present during 
inflation [56]. Because the asymmetry is not observed on small scales, both A and /nl are 
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required to be strongly scale dependent, which if parametrised by a power law would need 
to scale like [57] 

A oc /nl oc k nf NL , n/ NL < -0.56. (3.2) 


Such a large scale dependence normally invalidates the assumptions going into the usual 
calculation of nf NL , Eq. (1.2), if applied over more than about 2-4 e-foldings, making this 
parametrisation inadequate for studying the power asymmetry. For example, one expects 
dn/ NL /dlnk ~ n / NL when nf NL is larger than the slow-roll parameters [14, 15, 24], in which 
case n/ NL does not remain constant over much more than one e-folding when n/ NL ~ 1 . 

There is a strong tension between the requirement for a sufficiently large non-Gaussianity 
to explain the asymmetry and the stringent Planck constraints on both the bispectrum and 
the low l multipoles [58, 59, 49, 51, 60, 61], although a joint analysis of scale dependent /nl 
and scale dependent A should really be made to properly quantify the tension. Ignoring 
the scale dependence, the tension is quantified by the following equation (see e.g. [49] for 
the derivation), which relates A and /nl with the dipolar-induced quadrupole asymmetry, 

Cf } : 


c. 


(A) 


/nl 


4.2 x 10 


-li 


10 


= 6.6 


A 

007 


(3.3) 


where we have inserted the observed value of the squared quadrupole temperature anisotropy. 
While this does not have to correspond to the quadrupole generated by A, it can only be 
smaller if there is a finely tuned cancellation between the intrinsic quadrupole and this 
additional contribution. 

We are therefore motivated to study scenarios in which /nl is large (order of 100 [49]) 
on large scales l < 64 but reduces by an order of magnitude on small scales l > 700, which 
may be easier to arrange by making /nl cross zero on small scales. This corresponds to 
roughly 7 y NL ~ 5. The Planck forecast to be able to measure an error bar on n/ NL of [21] 


%L 


0 . 1 - 


50 

/nl 


(3.4) 


suggests 7 / nl ~ 5 may be detectable by Planck. This forecast was made assuming a fiducial 
value of /nl = 50 and constant n/ NL however, and so the analysis needs to be remade with 
the new parametrisation and new observational constraints. In order to do this, the full 
shape dependence of the bispectrum is required. 


We denote the super horizon mode (really the variation of a over the scale of the 
horizon) by A a, which is larger than a “normal” mode by an enhancement factor E, 
i.e. A a = EH^/{2 -k). This long wavelength mode will act as a shift in the classical back¬ 
ground value of the curvaton when the modes inside our horizon exit, meaning that opposite 
sides of today’s horizon will be sensitive to N a for different values of This concept was 
described in the context of inhomogeneous non-Gaussianity in [62, 63, 64, 65, 66 ]. Using 
a 5N expansion, and assuming a single-source model for the perturbations (meaning that 
only one field, such as a curvaton generates the primordial perturbation), we have up to 
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second order in Act 


2 A 




N a (a + A ff ) 2 - IV 2 


iV 2 


N a ((T + Act) — N a + N aa /S.(j + —N crcrcr A(T~ , 
= ^/nlAVAct — 2/nl^a/^C — /nl-EIO -4 , 


(3.5) 

(3.6) 


where 


N„ = 


dN 
d a k 


v ' = Ni 


_ 5 

/nl _ 6 m 


(3.7) 


Requiring |/nl| < 10 2 , which is reasonable since we should consider the constraint on 
/nl from large scales alone, we require E > 10, i.e. a super horizon mode whose amplitude 
is about 10 times larger than the typical sub horizon scale mode. The variation of the 
power spectrum due to long wavelength modes is measured by the squeezed limit of the 
bispectrum, which is closely related to local /nl- From (3.6) we can see that /nl ~ 10 3 
would naturally lead to A ~ 0.1, but the observational upper bound on /nl means that 
we instead need an unexpectedly large super horizon perturbation, E > 10, to generate a 
sufficiently large modulation of the power spectrum. 


3.1 The large trispectrum and quadrupolar power asymmetry 


When finding a regime with large and scale dependent /nl, we must also check whether 
<7 nl is within observational bounds. From Eq. (2.6), we have ^nl = 25g osc /(24r^ cc ) to 
leading order in l/rd ec - We are interested in regimes where |/ osc | <C 1, see (Eq. 2.6), 
but these do not correspond to regions where (?nl is accidentally suppressed so we should 
expect g nl 3> /nl, see e.g. Fig. 3 of [24]. It is straightforward to extend the analytic 
approximation of Sec. 2.2 to $nl- To first order in In (k/ko), we find in the limit r^ec 1 
and Sfc <C 1: 


S'nl(^) 


25 


2pL 


na r 


n (n 
~2~{n 


3) 

1) 


r] ao In (k/k 0 ) 


(3.8) 


Here, as before, ko corresponds to the scale at which /nl = —5/3, and a n = (n — 2)/(n — 1). 
For rj ao < 0.1, we see that over the observable range of scales $nl is approximately scale 
invariant: $nl ~ ~ n / r dec- Using Eq. (2.18) to substitute for rd ec we find 


Is'nl] ~ 4 x 10 5 



10 -2 


Vao 


2 


(3.9) 


This almost saturates the current observational bounds on ^nl, which are roughly |$nl| 

5 x 10° at 2 -ct [67, 68, 32, 69, 70] and recently tightened by more than a factor of two 
with the 2015 Planck results [71]. The problem can be alleviated for larger values of rj^o, 
however this will typically cause tension with the spectral index measured at the pivot 
scale ko, which for the self -interacting curvaton scenario is given by 


n°Q - 1 = — 2e^ 0 + 277^ 


(3.10) 
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Furthermore, the term quadratic in In A; becomes relevant for large rj ao (see Eq. (2.17)), at 
which point Eq. (3.9) breaks down. It is possible to generate a strongly scale-dependent 
/nl yet keep r/ CTo 10 2 if one considers > 1. In this regime, the spectral index can 
be kept compatible with observational constraints, however our numerical investigations 
confirm that an unacceptably large c/nl is still generated. As an example, we show in 
Fig. 4 a realisation of the self-interacting curvaton model with s/- > 1 in which /nl varies 
linearly from 30 on the largest scales to —10 on the smallest scales, crossing zero in the 
middle of the range of scales where the asymmetry is not observed. Without a dedicated 
analysis it is unknown whether such an /nl is compatible with the Planck constraint on 
constant |/nl| 1$ 10, but the figure of /nl as function of ( max from Planckl3 suggests it may 
be [72]. The scale dependence given in this example corresponds to 7 / NL ~ 5, which we 
can accurately treat as being close to constant over the observed range of scales provided 
that max[e^, 7 CT ] -C 0.1, see Eq. (1.4). The observed value of the spectral index implies 


l (multipole) 



Figure 4: The scale dependence of /nl oc A (blue line) and <?nl (red line) for n = 6, so — 3.85, 
r/ aa = 10~ 2 and rd ec = 0.05. The left-hand y -axis corresponds to the values of /nl and 10 _5 <7 nl, and 
the right-hand y-axis to the corresponding values of the asymmetry parameter A = 10^ 4 /nl-E-/ 2, 
for E — 42. Whilst the scale dependence of /nl is able to explain the CMB power asymmetry, 
the large value of ^nl — —2 x 10 6 is ruled out. This large value of $nl induces an unacceptably 
large quadrupolar asymmetry of the power spectrum: From Eq. (3.12) we find B ~ —14, which 
is almost four order of magnitude larger than the constraint B = 0.002 ± 0.016 of Ref. [73]. In 
this figure, /nl and $nl were computed numerically, with V(<f>) = m|/ 2 /2, = 16. To convert 
between wave number k and multipole l we have used the approximate relation l « x\ ss k, where 
^lss ~ 1-4 x 10 4 Mpc is the distance to the surface of last scattering. The pivot scale has been chosen 
to match the Planck pivot scale kg = 0.05 Mpc^ 1 . 


- 13 - 


A(k) = 10- 4 / NL (fc)F/2 







that this inequality should be satisfied and hence /nl will be linear in In k. Notice that 
a constant 7 / NL is incompatible with a scale dependent /nl and constant n/ NL , given the 
relation 7 y NL ~ /nlhj nl , showing that our new parametrisation is much more accurate in 
cases where ti y NL ~ 1- We also remark that there may be a large running of the spectral 
index due to the scale dependence of rj a (a correlation between this and large ^nl has been 
recently explored in [74, 75]), but this can be alleviated by making the contribution of rj a to 
the spectral index in (3.10) subdominant. We have numerically checked that the running 
is much smaller than the observational bounds in all of our explicit examples. 

A large <?nl is associated with a large induced quadrupolar asymmetry of the power 
spectrum. We parametrise this by B , 


P( = Piso + 2Ah.p + B (n.p) 2 ) , 


and is given by 


B = 


54 

xv3NL+ 



(JV ct Act ) 2 ~ 2g NL E 2 V ( 


4fl , NL-E' 2 * * 10 9 , 


(3.11) 


(3.12) 


where 


54 

25 


Isnl| 


N„ 


m 



(3.13) 


No direct observational constraint on B has been made, but an essentially equivalent 
quantity was constrained in the Planck2013 analysis of non-Gaussianity, see Fig. 17 of 
[72], where the data point with L = 1 corresponds to a rescaled version of A and the data 
point with L = 2 corresponds to B (where both parameters are assumed to be constant 
in the data analysis) 1 . This analysis did not require the preferred direction associated to 
A and B to be that same, unlike the prediction of our model but doing so should only 
make the constraint tighter. As an order of magnitude estimate, we require \B\ < 10 -2 
but none of our conclusions are sensitive to the exact value 2 . We note that there could 
be an independent contribution to B if there is more than one large, super horizon scale 
mode but one would have to tune both the amplitude and direction of the second mode 
in order to arrange a cancellation to the term we calculated here. For the example given 
in Fig. 4, where E rsj 40 and | 3 nl| rsj 10 6 * * we find \B\ ~ 14, which is about four orders 


1 We thank Zac Kenton and Antony Lewis for a discussion about this point. A recent paper by Kenton 
et al considers the possibility that asymmetry is generated by the next (third) order term which would 
appear in (3.11) [76]. The notation in that paper is related to ours by A = Ai and B = 2 A 2 . 

2 A comparable constraint was found on the related parameter g * defined in Fourier space 

P = Piso (1 + p*(n.p) 2 ) . (3.14) 

on which Kim and Komatsu [73] found the constraint g * = 0.002 ± 0.016, (see also [77] for the comparable 

Planck 2015 results). Their analysis did not require the direction of the quadrupolar power asymmetry to 

correspond to direction of the dipolar power asymmetry. Despite being defined in Fourier space, a preferred 

direction in Fourier space will generate a preferred direction on the CMB due to our view of it being a 

projection on a sphere. 
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of magnitude too large. While we cannot conclude that this is a problem for all models, 
and there are many free parameters to play with, we caution that parameter regions with 
large and strongly scale dependent /nl nray be associated with too large values of $nl 
and B, and their values should be checked whenever constructing an explicit model. A 
relation between large $nl and regions of scale dependent /nl was also spotted in [74] in 
the context of models with an observable running of the spectral index, but as far as we 
are aware the association with large B has not been made before. 

Another way to see why B is typically too large is remove the E dependence from 
Eq. (3.12) using Eq. (3.6), 


B ~ 1(T 2 


ffNL 

/nl(A0 2 



2 


(3.15) 


where care should be taken to evaluate A and /nl at the same pivot scale, since these terms 
are strongly scale dependent but B is not. Even for large self-interaction strengths, where 
no analytic solution exists, we see from Eq. (2.6) that regions where /nl crosses zero quickly 
(i.e. where 7 j- nl is large) will have r^ec 1 and |/ osc | 1- Such parameter regions do not 
correspond to an accidentally suppressed g 0sc , so we generically expect \B\ 2> 0.01 in the 
self-interacting curvaton scenario. It would be interesting to see if this conclusion applies 
to other models generating local non-Gaussianity. We remark that Eq. (3.15) applies to 
all single-source models (those in which perturbations from only one field generate the 
primordial curvature perturbation) and is not restricted to the self-interacting curvaton. 

Our parametrisation of /nl growing logarithmically is related to models in which a 
large “loop” generates the dominant contribution to /nl, which arises in cases where the 
tree level term to the bispectrum is suppressed, and such models also tend to have a large 
trispectrum [78, 79, 80]. However in such cases the role of the cut-off momenta is not fully 
understood [81], while our results are manifestly invariant of the pivot scale (which replaces 
the cut-off scale). See also the discussion in Sec. 3.2 of [24], and for a related perspective 
on the fine tuning required to have |5nl| |/nl| see [82], 

A curvaton potential which allows r] a < 0 may be preferred due to the observed value 
of the spectral index, and there has been interest in the axionic curvaton which does also 
generate a scale-dependent /nl [15]. However in this model the scale dependence is closely 
related to the running of the spectral index, which is tightly constrained by data to be 
small [83]. Hence the axionic curvaton cannot have a sufficiently strongly scale-dependent 
non-Gaussianity to explain the scale dependent power asymmetry. 

Overall, we conclude that there is no model of the dipolar power asymmetry which has 
been demonstrated to meet all observational requirements. The requirement for a strong 
scale dependence to this modulation means that the usual formalism of constant n/ NL will 
not be valid over the observed range of scales, while the formalism we introduced in this 
paper of yj NL = dlnfc/d/NL typically remains constant over the observed eight e-foldings, 
even when |7/ NL | 1. Regardless of the parametrisation used, we have shown that there is 

a high danger that $nl and/or B will be too large to match observations and their values 
should be checked when constructing models. 
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4. Conclusions 


In this paper we have introduced a new parametrisation for the scale dependence of /nl, 
which unlike previous parametrisations remains valid when /nl changes sign, and when 
the scale dependence is large. The parametrisation is simply given by 



(4.1) 


where 


d/NL 


(4.2) 


7/nl dlnfc' 


We expect 7 / NL to be approximately constant, and hence any higher order terms in Eq. (4.1) 
to be negligible, over a range which spans \n.(k/ko) inversely proportional to the slow-roll 
parameters, which is a much larger range than the seven e-foldings that we are able to 
constrain with CMB and LSS data. 

Given the tight observational constraints on constant /nl, our parametrisation is of 
most interest in the case /$ L ~ 0, and hence where there is a significant chance of /nl 
crossing zero within the observable range of scales. Our formalism is only valid for an equi¬ 
lateral configuration of the bispectrum in Fourier space, where the derivative depends on 
the single characteristic length scale. In order to observationally constrain our parametri¬ 
sation, it is necessary to be able to calculate the full bispectral shape. We leave this for 
future work. 

We have demonstrated the applicability of our formalism to the self-interacting curva- 
ton scenario, a model in which even a small correction to the mass term from self interac¬ 
tions can lead to strongly scale-dependent /nl • We have found the first analytical solution 
which remains valid if the curvaton evolves during inflation and the inflaton reheating, 
which despite being restricted to the limit of small self-interactions, is valid within the 
sufficiently non-linear regime to accurately describe /nl where it crosses zero. As a by 
product, we find that such regimes of enhanced scale dependence and /nl crossing zero do 
exist even in the limit of a quartic self-interaction term, which had not been previously 
found due to past analyses neglecting the evolution of the curvaton before the universe 
becomes radiation dominated after the inflaton has reheated. For arbitrary values of the 
self-interacting strength, we have found (4.1) to be a good match to numerical solutions, 
which also match the analytical solutions in the appropriate limit. 

We have applied our formalism to the observed power asymmetry of the primordial 
power spectrum. The asymmetry is required to be strongly scale dependent, since it is 
detected to be about a 10% effect on large scales l < 62, but must be at least an order 
of magnitude smaller on small scales, l > 600. This requires a similarly large decrease in 
/nl, something which our parametrisation can cope with but a commonly used previous 
one of constant n/ NL = din |/NL|/dln/c can not. The reason is that the previous formalism 
required |n/ NL ln(fc/fco)| 1$ 1, while observations require n/ NL < —0.5. 

While we can easily fine-tune curvaton models to have a sufficiently strong scale- 
dependent /nL) we find there are two generic new challenges which such a model must face. 
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In parameter regions where /nl is close to crossing zero and strongly scale dependent, we 
find that both <?nl and the quadrupolar power asymmetry are generically large and close 
to scale invariant. For the many examples we have tested, B and typically also (?nl are too 
large to agree with observations. In addition, there is the standard problem of explaining 
the origin of a very large amplitude super horizon scale perturbation which the power 
asymmetry requires. 
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A. General Expression for 


Here we provide the general dN expression for the bispectrum running parameter 7 / NL . A 
sum over /, J, K ,... is implicitly assumed. We start with the full expression for /nl: 


_ 5 NjjNjNj 

/NL 6 {N K N K y ’ 


(A.l) 


where Nj = dN/dipi k , etc and subscript k denotes evaluation of a quantity at the epoch 
of Hubble exit of mode k. Assuming slow-roll at horizon-crossing for all fields ipi we can 
write 

d _ <Pi k d 


d In k H k dtp Ik ’ 


(A.2) 


where the sum is carried out over all fields ipj. This enables us differentiate Eq. (A.l) with 
respect to k: 


7/nl — 


d/NL 


k=ko 
cO / 0 


din A: 

= — 2/nl(^ — 1 + 2eo) + 
where the spectral index reads: 


6 Ho 


Nij K N I Nj(p Ko 2 NjjN ik Nj^ Ko 


C n l n l y 


(n l n l y 


,(A.3) 


„o I _ o r , 2 VJoNijN! 

„ c _ i __ 2eo + 


(A.4) 
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and a superscript/subscript 0 denotes evaluation at the pivot scale k^. Using ^ = —Hp. 
and the slow-roll field equations, one can show that 

NjW Io = Wo , (A.5) 

NjjW Iq = W Jo - NjW IJo , (A. 6 ) 

N IJK W Io = W JKo - N u W IKo - N IK W IJo - NjWjjk 0 , (A.7) 


where Eqs. (A. 6 ) and (A.7) are derived by differentiating Eq. (A.5) with respect to c pj 0 . 
These expressions allow us replace third derivatives of N with second derivatives yielding 
the final result 


7/nl 


2/NL( n c — 1 + 2eo) 


5 

6 {N L N L f 


2 - Ar]iK 0 NijNjN K 


Vtj 0 NrNj - £ijk 0 NiNjN k , 


(A. 8 ) 


where r//j = Wu/W and f jjk = Wijk/W. In the limit where the primordial perturba¬ 
tions are generated by a curvaton field, a, (the single source limit) the expression for y^ NL 
simplifies to: 

wO c: wO 

^ = inf sk " w °' 3r ‘ /2 (single source) ’ (A ' 9) 

where r — 8/N% is the tensor-to-scalar ratio. Although r must be small, W^ aa /H q is not a 
slow-roll quantity and it may be much larger than unity in models with a non-quadratic 
potential. This result is obtained by carefully taking the limit 


f^h N J 
\ e N 

y c cr o 1 v cr 


< l, 


VJJo 

Vcrcro 



(J 


(A-10) 


in Eq. (A. 8 ). One can easily verify that inserting N a = (2/3 )r<i ec (o' osc /a osc ) with a osc given 
by Eq. (2.12) into the above equation, one obtains the result for 7 f NL given in Eq. (2.18), 
as expected. 
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